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We propose a method to extend into the bulk asymptotically flat static spherically sym- 
metric brane- world metrics. We employ the multipole (1/r) expansion in order to allow 
exact integration of the relevant equations along the (fifth) extra coordinate and make 
contact with the parameterized post-Newtonian formalism. We apply our method to 
three families of solutions previously appeared as candidates of black holes in the brane 
world and show that the shape of the horizon is very likely a flat "pancake" for astro- 
physical sources. 
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1. Introduction 

In recent years models with extra dimensions have been proposed in which the Stan- 
dard Model fields are confined to a four-dimensional world viewed as a (infinitely 
thin) hypersurface (the brane) embedded in the higher-dimensional space-time (the 
bulk) where (only) gravity can propagate. Of particular interest is the Randall- 
Sundrum (RS) model with one infinitely extended extra dimension "warped" by 
a non- vanishing bulk cosmological constant A 032 In this context, it is natural to 
study solutions corresponding to compact sources on the brane, such as stars and 
black holes. This task has proven to be rather complicated and there is little hope 
to obtain analytic solutions such as those found with one dimension less ^1. The 
present literature does in fact provide solutions on the brane 4 5 SI perturbativc 
results over the RS background EEI and numerical treatments EES 

In this letter we investigate how to extend into the bulk asymptotically flat static 
spherically symmetric solutions on the brane. We hence set the brane cosmological 
constant to zero by fine tuning A to the brane tension, which we denote by 3 a (in 
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units with 8 it G( $) = 1, where G/g) is the five-dimensional gravitational constant), 
i.e. A = —3 u 2 /2 1 1 ! 11 ! We then note that the bulk metric in five dimensions must 
satisfy {fi, v = 0, . . . , 4) 

WR„v = -Ag MV . (1.1) 

On projecting the above equations on the brane and introducing Gaussian normal 
coordinates x l (i = 0, . . . , 3) and z (z = on the brane), one obtains the constraints 



= Mr 

z=0 



= 0, (1.2) 

z=0 



where ^R is the four-dimensional Ricci scalar and use has been made of the neces- 
sary junction equations One can view Eqs. (|1.2f) as the analogs of the momentum 
and Hamiltonian constraints in the Arnowitt-Dcser-Misner (ADM) decomposition 
of the metric and their role is to select out admissible field configurations along a 
given hypersurface of constant z. Such field configurations will then be "propagated" 
off-brane by the remaining Einstein equations, namely 

' ""//,, A'/,, . (1.3) 
The above "Hamiltonian" constraint is a weaker requirement than the purely four- 
dimensional vacuum equations Rij = and is equivalent to Rij — E^j , where Eij is 
(proportional to) the (traceless) projection of the five-dimensional Weyl tensor on 
the brane 

ED 

In the next Section we outline a general procedure to solve the bulk equations 
for a given static spherically symmetric brane metric, which we shall then apply to 
several cases in Section [21 In Section 0] we finally summerize and comment on our 
results. 

2. The procedure 

We consider five-dimensional metrics of the form 

ds 2 = -N(r, z) dt 2 + A{r, z) dr 2 + R 2 {r, z) dQ 2 + dz 2 (2.1) 

where dQ 2 = d6 2 + sin 2 9 d<p 2 and N, A and R are functions to be determined. For 
such metrics the momentum constraint is identically solved and the "Hamiltonian" 
constraint reads 



4 



where the subscript B is to remind that all functions are evaluated on the brane at 
z = 0, ' = d/dr and we set Rb = r thanks to the spherical symmetry 
Our strategy to construct the bulk metric is then made of three steps: 

(i) choose a metric of the form p.lfl whose projection on the brane, 

ds\l= - -N B (r) dt 2 + A B (r) dr 2 + r 2 dtf , (2.3) 

solves the constraint Eq. 1)2. 2fl : 



I 



I 
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(ii) expand such a metric in powers of 1/r (four-dimensional multipole expansion) 
to order n, i.e. introduce the quantities 

N n (r, z)} "if 

A n (r,z)\=J2^\^(z) , (2.4) 
Rl{r,z)) fc=o [r 2 c k (z) 

where n/c(0), a&(0) and Cfe(0) reproduce the particular solution chosen at step 1 
(to order n): 

n fc (0) ^ (N B (r) 



k— o ^ r Cfe(0) 



(iii) substitute the sum l|2.4l) into Eq. 1)1. 3|) and integrate exactly in z the (three) 
equations thus obtained for the functions n n (z), a n (z) and c n (z). 

This procedure turns out to be particularly convenient for the problem at hand 
because it converts Einstein equations l|1.3l) into three sets of second order ordinary 
differential equations (in the variable z) of the form 



dz 2 



- o l f n = F k<n , (2.6) 



where /„ is any of the functions n n (z), a n (z) and c n (z) (n > 1), and Fk< n {z) a 
functional of the lower order terms /fe< n 's and their first and second derivatives. 
The relevant boundary conditions for Eq. 1)2. 6J) are given by the requirement 1)2.5(1 
for /„(0) and the junction conditions which imply d/„/dz(0) = — af n (0). We 
shall not display the (cumbersome) Fk< n for simplicity, however they turn out 
to be such that the Cauchy problem thus defined admits analytic solutions. The 
entire procedure can then be executed automatically with the aid of an algebraic 
manipulator to determine the functions /„ recursively, from the lowest order up a , 
with the only limitation of the power and memory of the available computer. 

Let us now comment on a few more points regarding the above procedure. First 
of all, we wish to stress that there is a large freedom in the choice of the metric on the 
brane. In particular, the coefficients rtfe(0) and afe(0) can be chosen at will, except 
for the algebraic constraints following from Eq. 1)2. 2J) . Since such coefficients are 
related to the shape of the source, this input represents the physical content of the 
model. A second, related point is the convergence of the series expansion 1)2.4)1 . It is 
in general difficult to pinpoint one parameter (among the many possible coefficients 
of the multipole expansion) whose "smallness" guarantees that orders higher than 

a For n = one has a system of three coupled second order ordinary differential equations for the 
fo's. The corresponding Cauchy problem is solved by the usual warp factor, /o = exp(— a z), which 
is unique as follows from the usual theorems of uniqueness. We further note that for r — ► oo our 
bulk solutions reproduce the RS space-time Q as one expects for an asymptotically flat brane 
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n be negligible. Because of this, we should consider our results as reliable for those 
values of r and z such that 

yn-\-l 

for given values of the parameters rifc(O) and cifc(O). In general, for a given z, such 
a condition will be satisfied for sufficiently large r. 



E 

k=0 



fk(z) 



(2.7) 



3. Application to brane-world black holes 

As examples of brane metrics, we have considered the solutions given in Refs. I^ISEl 
which can be expressed in terms of the ADM mass M and the post-Newtonian 
parameter (PNP) measured on the brane. The case wi th 77 = (exact 

Schwarzschild on the brane) is the well known black string (BS)^l which extends 
all a long th e extra dimension. The BS is known to suffer of serious stability prob- 



lems' 



e.g. the Kretschmann scalar, 



K 2 = WRpuXp ^R^ xp = a 4 + 48 M 2 



(3.1) 



diverges on the AdS horizon (z — » 00). One is therefore led to conclude that black 
holes on the brane must depart from Schwarzschild and have 77 ^ 0. As was suggested 
in El the interesting cases are those with 77 < 0, since 77 > implies some sort of 
anti-gravity effects (see later for further comments). 

Short distance tests of Newto nian gravity yield the bound cr _1 < lmm^and 
I77I < 10~ 3 from solar system tests ^4 Since we want to study astrophysical sources 
of solar mass size, in the following we shall often refer to the typical values 



M = 10 a 



7 „-l 



1 km 



V 



-10" 



(3.2) 



In such a range (Ma ^> 1 an d M-C 1) one finds a qualitatively identical behavior 
for all brane metrics in Refs. I ^ I Q I so we shall just give the results for case I of 
Ref.©(see also El), that is 



Ab = 1 



An = 



2M 
r 



3M 
~2r 



(3.3) 



2M 



3M 
~2r 



9 



V 



where r = = 2 M is the event horizon and the remaining (non- vanishing) PNPs 
are (3 = j = 1 + ^rj. 

We applied the above procedure to the brane metric l|3.3|) and were able to solve 



I 



I 
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the corresponding Eqs. (|2.6|l up to n — 19. For brevity, we just display a few terms: 



N, = e 



A, 



2M r]M 2 (1 - e° 



r 3r 4 

2M , 
1 + — 3 + 77 
3 r 



M 
1 + — 

r 



AM 2 ( 77 77 2 



7- 



8M 3 2t]M(l 



3r 3 



(3.4) 



i + 



rjM (1 - e CTZ ) 
3 r 3 a 2 



2Af 7M 2 



2r 2 



It is important to note the appearance of positive exponentials in the metric func- 
tions. Such terms (which also show up at higher orders) are non-perturbativc in z, 
which makes the expansion in 1/r preferable (or at least complementary) to the 
expansion for small z. 

For 77 < one generally finds that, for every n > 3 and r > 0, there exists a 
corresponding value of z [say z^ xls (r)) such that i? 2 (r, zf* ls ) — 0. Since An R 2 is the 
proper area of the sphere t = r = z — constant, this seems to indicate that the axis of 
cylindrical symmetry is given, in our Gaussian reference frame, by a line z = z axls (r) 
(which should be exactly obtained in the limit n — > 00). Since i? 2 is determined 
up to corrections of order n + 1, which in general do not vanish for z — z^ xls (r), 
one cannot consider this as a mathematically rigorous proof [the condition l|2.7|l 
obviously fails for z = z^ xls ]. However, we found that in the physically interesting 
range of the parameter space [M, 77) the 1 jr expansion yields rather stable values 
of z^ xls (r) in a wide span of n. The stability improves for larger values of r [as could 
be inferred from Ij2.7|l ] and becomes very satisfactory for r > (see Fig.^). From 
Eq. I|3.4|l and Ma > 1, r > r h , one finds 



1 , /3a 2 r 3 
- — In — 

2(7 \-T]M 



(3.5) 



which numerically agrees fairly well with zfg IS (r). Further, we have checked that 
R n (ri,z) > R n (r2,z) for n > r-2 and < z < z^ xls (r2), so that lines of constant 
r do not cross and the Gaussian normal form 12.1fl of the metric is preserved for 
77 < in the bulk within our approximation (for an example see Fig.|21). 

It is interesting to compare our solutions for 77 < with the BS . In particular, 
one would like to see the shape of the horizon in the bulk, knowing that for the BS 
it does not close but extends all the way to the AdS horizon. First we note that, if 
the horizon closes in the bulk, then it must cross the axis of cylindrical symmetry 
at a point (the "tip") of finite coordinates (r tlp , z tlp ) where = R 2 = 0. Of course, 
we just have such equations explicitly at order n, 



N n (r^,z^) = Rl(r« p ,z^) = 



(3.6) 
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For large values of n, one can solve Eqs. I|3.6I) numerically and find the "tip" . More 
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11 12 13 14 15 16 



z 



Fig. 2. The function R 2 evaluated to order n = 19 for r = M, r\y = 2M (dashed line) and 10 M 
with M = 10 6 a' 1 and r) = -1CT 4 . 




Fig. 3. Qualitative picture of the bulk structure. 

in detail, for n = 19 one has 

N 19 (r h + e : z) > and N lg (r h - e, z) < , (3.7) 

f° r £ ^ (3/100) rh and < z < 0.99 zj l g ls (r). Thus, to a very good approximation, 
N n (r^,z) ~ for < z < 2^ xls (rh) when Mc 3> 1 and |?y| <C 1 and negative. 
A good parameterization for the horizon is thus given by r ~ and < z < 
■ 2 ^ xls ( r h) — z* lp (see Fig.|3J|. For the typical parameters l|3.2ll we have r^f ~ 0.97 
and zf g p ~ 20.8 cr -1 which is very close to zfg ls (rh) = 20.9 a^ 1 . This all strongly 
suggests that the horizon does close in the bulk, as previously obtained by numerical 
analysis 1 ^ 1 -^ 1 (for a comparison with the BS see Fig.^J). 

One can now get an estimate of how flattened the horizon is towards the brane 
by comparing the proper length of a circle on the horizon which lies entirely on 
the brane, Cm = 2nrh ~ 1.3 • 10 s c -1 , with the length of an analogous curve 
perpendicular to the brane, C± ~ 4z' g p ~ 84 cr _1 . Since their ratio is huge, one can 
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Fig. 4. The function R 2 /R 2 jS evaluated to order n 
10 M and 100 M with M = lO 7 ^- 1 and r? = -10" 4 . 



19 for r = M, r h = 2M (dashed line), 



in fact speak of a "pancake" horizon as was suggested, e.g., in Ref. 

It is interesting to note that for n = 5, one can still solve Eqs. (|3.6(l analytically 
and finds 



tip 



1 



= - In 
a 



1 



a r n a 

V-v 



(3.8) 



where a ~ 0.5 and we used Afcr > 1 in the final expression. This yields Zg lp = 
20.7 <7 -1 for the parameters 13. 2|) . in excellent agreement with the numerical value 
obtained at order 19. In light of this stability, one can therefore approximate the 
dependence of the exact z tlp on the black hole ADM mass M from Eq. (|3.8I) and 
obtains that the area of the (bulk) horizon is approximately equal to the four- 
dimensional (brane) expression b , 



(5 U = 4tt / R z (r h ,z)dz 



16 7T 



(2M) 2 



(3.9) 



where we again used Ma 3> 1. Eqs. I|3.8II and 1)3. 9|) again supports the idea of a 
"pancake" shape for the horizon (see also^ for the logarithmic dependence of z tlp 
on M). 

Drawing upon the above picture, in particular the crossing of lines of constant r 
with the axis of cylindrical symmetry at finite z, one can infer that the space-times 
we obtain do not suffer of one of the instabilities of the BS, namely the diverging 
Kretschmann scalar E3. i n f ac t ; K 2 is still an increasing function of z along lines of 



The fundamental (possibly TeV scale) five-dimensional gravitational coupling G( 5 ) ~ Giv/c, 
where G^r is the four-dimensional Newton constantly Thus, from 13.91 one has ( 5 )«4/G( 5 ) ~ 
AI 2 /G N ~ < 4 U/Giv. 



I 



I 
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constant r, but one has 

K 2 -^ = a^+o(^<l^ + o(^) , (3.10) 

where we used Eq. (|3.5|l to maximize K 2 uniformly in z. The coefficients a and b 
depend on the parameters of the multipole expansion, correctly vanish in pure RS 
and b = 48 M 2 for r\ = [the BS, cfr. Eq. I|3.1|) ]. The remaining problem of stability 
under (linear) perturbations is a difficult one to tackle and we do not attempt 
at it here. 

As we mentioned previously, the cases with r\ > show a very different qual- 
itative behavior. One in fact finds that R^(r, z) is generically a (monotonically) 
increasing function of z for all (su fficie ntly large) values of r, as one would indeed 
expect on a negative tension brane^^. However, for any r±, r-i > there now exists 
z n( r i> r 2) such that i?^(ri,z*) = i?^(r 2 , z*), i.e. space-like geodesies of constant r 
display caustics and the Gaussian coordinates (r, z) do not cover the whole bulk EH 

4. Conclusions and outlook 

We have explained in some detail a method to extend into the bulk a given asymp- 
totically flat static spherically symmetric metric on the brane which is based o n the 
multipole (1/r) expansion. The application of our method to candidate metrics HEED 
for astrophysical sources led us to conclude that black hole horizon closes in the bulk 
and indeed has the shape of a "pancake" . Our solutions depend on three parameters 
(a, M and 77), although one could argue that only two of them are independent. For 
instance, in order to recover the four-dimensional Schwarzschild metric when the 
extra dimension shrinks to zero size (a -> 00 in the brane equations El) , one can 
guess 1 77 1 ~ l/M a. Reducing the number of dimensions is however a singular limit 
of the five-dimensional metric, so it is critical to obtain precise relations among the 
parameters by this procedure. It is also difficult to extract sensible results for tiny 
black holes (Ma < 1), e.g. from Eq. one now gets z" p ~ M and Eq. 

yields the relation for five-dimensional Schwarzschild black holes ~ M 3 LLU 

Our expansion however suggests that the horizon departs significantly from the line 
r = ?"h when M a -C 1 and the above estimate is very rough. The dependence of the 
horizon area on the ADM mass is crucial to study the Hawking evaporation and we 
hope to return to it in the future. 
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